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PREFACE 

This  final  report  summarizes  a  research  program  performed  by  the 
McDonnell  Douglas  Research  Laboratories,  St.  Louis,  Missouri,  on  the  response 
of  ramjet  inlet  flows  to  periodic  downstream  disturbances.  The  disturbances, 
generated  by  a  mechanical  exciter  In  cold  flows,  simulate  pressure  fluctua¬ 
tions  that  can  be  exhibited  by  dump  combustor  Instabilities  in  tactical-size 
ramjets.  The  research  was  conducted  under  Contract  N00014-80-C-0481  for  the 
Office  of  Naval  Research.  The  performance  period  was  1  July  1980  to  31 
January  1983. 

The  principal  investigator  was  Dr.  Miklos  Sajben;  Dr.  Thomas  J.  Bo gar 
and  Mr.  Joseph  C.  Kroutil  were  coinvestigators.  The  program  manager  was 
Dr.  Richard  S.  Miller,  Office  of  Naval  Research. 

This  report  has  been  reviewed  and  approved. 


R.  J.  Hakkinen 
Director-Research 

McDonnell  Douglas  Research  laboratories 


D.  P.  Ames 

Staff  Vice  President 

McDonnell  Douglas  Research  Laboratories 
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NOMENCLATURE 


A  cross-sectional  area 

B,E  Fourier  coefficients 

f  frequency 

h  channel  height 

i,j  dummy  Indices 

M  Mach  number 

M  average  Mach  number  defined  by  Equation  (2) 
n  dummy  Index 

p  static  pressure 

P  rms  value  of  time-dependent  pressure;  normalized  by  rms  value  of  exit 
amplitude 

R  perfect  gas  constant 

t  time 

T  absolute  temperature 

u  x-component  of  velocity  vector 

x  streamwise  (horizontal)  coordinate;  x  ■  0  at  ramp  lip,  increasing 
downstream 

y  transverse  coordinate;  y  -  0  at  ramp  lip,  increasing  upward 

z  spanwlse  coordinate;  zero  in  plane  of  symmetry  of  model 

y  ratio  of  specific  heats 

v  pressure  ratio  (lnlet-total/exit-statlc) 

€ 

p  density 

Subscripts 
«  freestream 

o  cowl  lip  (geometric  throat) 

o  shock 

e  exit  station 

v  choke  point  at  exhaust;  approximately 

coincides  with  rotor  axis 

1,2, 3... 6  designating  orders  of  harmonics 

c  critical 

b  buzz 

f  narrow-band  filtered  at  excitation  frequency 


x-location 
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1.  INTRODUCTION 


The  liquid-fuel  ramjet,  combined  with  an  integral  solid  rocket  booster, 
is  a  candidate  for  high-speed  tactical  missile  propulsion  systems.  The  ex¬ 
perience  of  development  programs  indicates  that  such  systems  may  display  low- 
frequency  «  300  Hz),  large-amplitude  (±  15%)  pressure  oscillations  that  in¬ 
volve  both  the  combustor  and  the  inlet.1  If  the  initial  shock  system  of  the 
inlet  is  sufficiently  disturbed  by  these  fluctuations,  the  inlet  may  unstart, 
resulting  in  a  large  reduction  of  captured  airflow  and  thrust,  and  eventually 
causing  failure  of  the  mission. 

Pressure  oscillations  in  dump-combustor  ramjets  have  been  investigated  by 
Hall,1  Rogers,2,3  Clark,4*5  Schadow  et  al.,6  Crump  et  al.,7  and  Waugh  et 
al.®  One  of  the  widely  accepted  conclusions  of  these  studies  is  that  the 
inlet  is  a  significant  element  of  the  problem  and  that  the  dynamics  of  inlets 
require  detailed  investigation. 

The  present  contract  is  intended  to  generate  such  information  through  the 
experimental  study  of  a  seml-freejet  inlet  model  subjected  to  periodic,  con¬ 
trolled  disturbances  at  its  downstream  end,  simulating  combustor-pressure 
fluctuations.  The  disturbances  are  introduced  by  the  periodic  mechanical  mod¬ 
ulation  of  the  choked  exhaust  cross-sectional  area. 

The  flow  in  the  same  inlet  model  with  steady  exit  boundary  conditions 
formed  the  subject  of  a  separate,  MDRL-sponsored  investigation,  in  which  both 
time-mean  and  the  time-dependent  aspects  of  the  flows  were  explored.  Documen¬ 
tation  of  the  results  is  available®;  familiarity  with  this  reference  is  help¬ 
ful  in  understanding  the  present  report. 
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2.  OBJECTIVES 


The  objectives  of  this  program  were  contractually  stated  as  follows: 


Task  1  - 

Task  II  - 


Task  III  - 


Design  and  build  an  exciter  device  for  periodic  modulation  of 
the  exhaust  area  of  an  appropriate  diffuser  model. 

Obtain  still  and  motion  photographs  of  unsteady  flowfields  in 
a  supercritically  operated  diffuser  over  suitable  combinations 
of  shock  strength,  excitation  frequency,  and  excitation 
amplitude. 

Determine  unstart  boundaries  representing  combinations  of 
excitation  frequency  and  amplitude  that  cause  the  shock  to 
pass  through  the  inlet  throat  for  the  available  range  of 
terminal  shock  Mach  numbers. 


Task  IV 


Determine  phase-averaged  shock  displacement  histories  for 
selected  combinations  of  operating  parameters. 


3.  DESCRIPTION  OF  EXPERIMENT 


3.1  Inlet  Model 

The  original  test  concept  involved  a  dual-throat  model  configured  as  a 
small  supersonic  wind  tunnel.  References  to  a  diffuser  in  the  task,  statements 
reflect  these  initial  ideas.  However ,  it  soon  became  clear  that  the  principal 
objectives  would  be  better  served  by  an  inlet  model  in  a  semi-freejet  arrange¬ 
ment,  and  the  approach  was  changed  accordingly.  An  inlet  model  available  from 
an  MDC-sponsored  investigation  was  used.  The  model  (Figures  1,  2  and  3)  is  an 
external-compression*,  ramp-type  inlet  operated  in  a  semi-freejet  mode  in  the 
exhaust  of  a  supersonic,  nozzle.  The  ramp  angle  is  14.7°,  which  in  conjunction 
with  the  freestream  Mach  number  of  1.84  produces  an  attached  weak  oblique 
shock,  yielding  a  ramp-surface  Mach  number  of  1.3.  This  Mach  number  is  close 
to  that  found  in  typical  ramjet  inlets ,  although  it  may  be  produced  by  a 
combination  of  higher  freestream  Mach  numbers  and  higher  ramp  angles.  The 
oblique  shock  is  designed  to  clear  the  cowl  lip  by  0.7  mm,  such  that  the  cowl 
lip  is  also  exposed  to  the  post-shock  Mach  number  of  1.3. 

The  inside  surface  at  the  cowl  lip  is  parallel  to  the  flat  ramp  surface, 
and  the  cross-sectional  area  increases  monotonically  beyond  this  point  so  that 
a  geometric  throat  (dh/dx  -  0)  is  formed  at  the  cowl  lip.  The  height  of  the 
channel  at  this  location  (hQ  -  23.5  mm)  is  used  as  a  reference  length  through¬ 
out  the  report.  The  design  makes  no  allowance  for  boundary-layer  displacement 
effects  on  either  surface;  the  flow  contains  a  weak  oblique  shock  (Mach  wave) 
initiated  by  boundary-layer  growth  at  the  cowl  lip.  The  channel  continues  to 
diverge  (maximum  angle  is  6.9°)  and  forms  a  gradual  transition  to  a  constant- 
area  segment.  The  contours  are  simple  combinations  of  straight  lines  and  cir¬ 
cular  arcs;  no  attempt  was  made  to  optimize  any  performance  indicator  by 
tailoring  the  wall  shapes. 


*Mixed  compression  inlets  are  often  used  in  ramjet  technology.  They  are  more 
efficient  at  high  speeds,  but  must  have  variable  geometry  for  starting.  The 
introduction  of  such  requirements  into  the  present  study  would  have  diverted  a 
needlessly  large  fraction  of  resources  from  the  central  issue,  and  was 
therefore  avoided. 
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All  dimensions  in  cm  ^V^Jnlet Exciter 


•  Inlet  and  exciter  width:  17.78  (constant) 

•  Aspect  ratio  at  throat:  7.56 

at  exit:  2.80 

GP3I4I14-2 

Figure  1.  Inlet  model  including  the  supersonic  nozzle,  inlet,  and 
exciter  components. 
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Figure  2.  Inlet  model  installed  in  test  facility. 
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Figure  3.  Inlet  model  contour  shapes  and  locations  of  fast-response  pressure  transducers  indicated  by 
solid  triangles.  Origin  of  the  coordinate  system  is  at  the  ramp  lip. 


An  arbitrarily  chosen  location  within  the  constant-area  terminal  segment 
is  designated  as  the  exit  station  (xg  *  21.91),  which  serves  as  the  downstream 
boundary  of  the  flowfield  of  interest  and  corresponds  to  a  station  just  up¬ 
stream  of  fuel  injection  in  real  ramjets. 

The  supersonic  flow  was  produced  by  a  two-dimensional,  convergent- 
divergent  supersonic  nozzle,  designed  by  the  method  of  characteristics  and 
corrected  for  boundary-layer  displacement  effects.  The  inlet  is  located  in 
the  core  flow  of  the  nozzle  exit;  the  top  and  bottom  nozzle-wall  boundary 
layers  are  removed  through  large  slots  above  and  below  the  ramp  and  cowl  lips, 
respectively.  Approximately  68Z  of  the  nozzle  flow  enters  the  inlet. 

The  nozzle  is  supplied  with  dry,  filtered  air  from  a  51-cm  diameter,  160- 
cm  long  plenum  chamber  immediately  upstream.  The  plenum  has  several  flow  con¬ 
trol  screens  and  a  smoothly  contoured  transition  to  the  nozzle  with  a  contrac¬ 
tion  ratio  of  22.4:1,  providing  highly  uniform,  low-turbulence  flow.  Plenum 
chamber  pressures  ranged  up  to  570  kPa  with  mass  flow  rates  through  the  nozzle 
of  6.7  kg/s. 
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The  nozzle  and  Inlet  form  a  single  structural  unit  with  common  sidewalls 
that  permit  glass  windows  to  overlap  the  gap  between  the  nozzle  and  the  Inlet 
This  design  assures  excellent  optical  access  to  a  region  in  which  Important 
processes  occur  during  start,  buzz,  and  unstart.  However,  the  lack  of  side- 
wall  boundary- layer  control  adversely  affects  two-dimensionality  and  a  large 
throat  aspect-ratio  (7.56:1)  was  chosen  to  minimize  sidewall  boundary-layer 
influence. 

3.2  Exciter 

Exit  area  modulation  was  accomplished  by  the  exciter  device  shown  in 
Figures  1-4.  The  principal  components  of  the  exciter  are  an  elliptical  cross 
section  rotor  with  axis  horizontal  and  normal  to  the  flow,  two  flaps  that  can 
be  tilted  around  the  indicated  hinge  points,  and  two  doors  formed  by  the 
hinged-end  segment  of  the  top  and  bottom  walls.  These  components  formed  a 
symmetric  arrangement  of  four  parallel  flow  passages,  each  shaped  as  a 
convergent-divergent  channel  and  choked  under  all  operating  conditions. 


Am  «  =  0-64*1 
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All  dimensions  in  cm 


Flap  letting  for 
mwdniam  ampHtnde 


v  Flap  setting  for 
zero  ampHtnde 


Fignre  4.  Mechanical  exciter  and  method  of  ampHtnde  control. 


The  total  cross-sectional  area  at  the  choke  point  (the  sum  of  the  four 
passage  areas,  referred  to  as  exhaust  area)  could  be  varied  according  to  the 
relation 

Ay  -  Ay  +  A'  cos(2*fext)  ,  (1) 

where  A  is  the  time-mean  value  and  A",  f  are  the  amplitude  and  frequency 
v  v  ex 

of  the  time-dependent  contribution,  respectively.  Equation  (1)  does  not 
describe  higher  harmonics,  which  were  present  with  negligible  amplitudes  only 
and  will  not  be  explicitly  considered. 

Variation  of  the  three  parameters  included  in  the  right  side  of  Equation 
(1)  was  continuous,  noninteracting  (or  at  most  weakly  interacting),  and  could 
be  accomplished  during  runs  without  shutting  down  the  experiment.  The  sym¬ 
metry  of  the  configuration  represented  an  improvement  over  an  asymmetric 
device  used  in  earlier,  similar  experiments.*®'** 

The  doors  and  flaps  were  actuated  remotely  by  appropriate  motors  and 
linkages  that  assured  symmetric  positioning  of  the  respective  component  pairs. 
The  rotor  was  driven  by  a  0.56  kW  variable-speed  motor  via  V-belt  pulleys. 

The  doors  fixed  the  height  of  the  two  outer  passages,  and  since  the  time- 
mean  area  of  all  passages  was  carefully  designed  to  be  independent  of  flap 
setting,  the  time-mean  area  (A^)  was  determined  solely  by  the  door  setting. 

The  doors  thus  performed  the  control  function  of  the  conventional  plug 
throttle  commonly  used  in  inlet  testing.  They  were  not,  however,  used  for 
flow  measurements  as  calibrated  plugs  usually  are,  for  two  reasons.  First, 
the  exciter  assembly  was  not  preceded  by  a  sufficiently  long  channel  segment 
to  create  a  flow  distribution  independent  of  far-upstream  events,  as  would  be 
required  for  a  reliable  flow-rate  vs.  open-area  calibration.  (Making  the 
channel  long  would  have  altered  the  natural  frequencies  of  the  system.) 
Secondly,  the  effective  time-mean  area  of  the  exciter  was  slightly  dependent 
on  excitation  amplitude  and  frequency,  introducing  additional  error  sources 
into  any  flow  measurement  attempt.  The  mean  exit  area  thus  did  not  have  a 
simple  connection  to  mass  flow  and  was  therefore  considered  inappropriate  for 
use  as  a  parameter  to  characterise  time-mean  flow  conditions. 
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The  rotor  nodulated  the  heights  of  the  two  loner  passages  at  twice  the 
shaft  frequency,  with  a  constant  amplitude  of  3*18  mm,  or  13. 5Z  of  the  geo- 
aetrlc  throat  height.  The  amplitude  of  the  nodulatlon  was  the  sane  for  both 
passages,  but  the  relative  phase  between  the  two  could  be  varied  by  using  the 
flaps.  In  the  flap  position  shown  in  Figure  4a,  the  two  areas  vary  in  phase; 
the  minima  and  maxima  occur  simultaneously  on  top  and  bottom.  If  the  flaps 
are  set  as  shown  in  Figure  4b,  then  the  area  variations  are  180°  out  of  phase, 
and  the  occurrence  of  a  minimum  on  one  side  is  associated  with  a  maximum  on 
the  other  side.  The  sum  of  the  two  inner  passage  areas  is  then  constant,  and 
the  modulation  amplitude  (Av')  is  zero.  Intermediate  flap  settings  produce 
intermediate  amplitudes . 

An  optical  encoder  mounted  on  the  rotor  shaft  provided  a  pulse  train  for 
measurement  of  the  rotor  speed  and  use  as  a  reference  in  ensemble-averaging 
various  signals.  When  needed,  steady  boundary  conditions  were  created  by 
locking  the  rotor  with  its  major  axis  horizontal,  such  that  the  inner  passage 
areas  were  maximized. 

After  an  initial  development  phase,  the  exciter  performed  reliably  and 
according  to  expectations. 

3.3  Instrumentation 

Fifteen  channels  of  steady  data  were  routinely  recorded.  Including  opera¬ 
tional  variables  such  as  pressures  in  the  plenum  chamber  and  at  various  loca¬ 
tions  along  the  top  wall  of  the  nozzle  and  inlet  model,  temperatures,  and  rms 
values  of  selected  surface  pressures. 

The  top  and  bottom  walls  of  the  model  contain  over  100  orifices  for  mea¬ 
suring  time-mean  wall  static-pressure  distributions.  As  shown  in  Figure  3  and 
Table  1,  10  ports  are  located  along  the  top  wall,  and  2  along  the  bottom  wall, 
to  accommodate  sensors  for  fluctuating  surface  pressures.  Two  similar  ports, 
one  each  in  the  cowl  and  ramp  slots,  were  used  to  detect  passage  of  the  ter¬ 
minal  shock  over  the  cowl  and  ramp  lips.  Miniature  straln-gage-type  trans¬ 
ducers  with  flat  frequency  response  to  90  kHz  were  used  to  measure  the  un¬ 
steady  pressures.  Surveys  of  the  unsteady  wall  pressures  were  amplified,  low- 
pass  filtered  at  10  kHz,  and  routed  either  for  recording  on  FM  tape  or  for 
further  analog  signal  conditioning. 


TABLE  I.  FAST-RESPONSE  TRANSDUCER  LOCATIONS. 


Transducer  symbol 
oa  Figure  3 

^0 

X 

Top-wall 

RS  (ramp  slot) 

0.473 

RI  (rump  inside) 

1.035 

T1  (top-1) 

2.519 

T2 

4.147 

T3 

5.761 

T4 

7.373 

T5 

8.995 

T6 

11.114 

T7 

13.829 

T8 

17.349 

ET  (exit  top) 

21.914 

Bottom-wall 

CS  (cowl  slot) 

1.862 

Cl  (cowl  inside) 

2.348 

EB  (exit  bottom) 

21.914 

cra-wM-i 


A  vertical  rake  containing  four  fast-response  pressure  transducers 
(1.6  mm  o.d.)  was  used  to  measure  the  time-dependent  total  pressures  at  the 
exit  station  midspan.  The  rake  could  be  stepped  in  the  vertical  direction  to 
three  fixed  locations,  adequately  resolving  the  vertical  variation  of  the 
unsteady  total  pressure  (12  points). 

High-speed  (5000  frames/s)  schlieren  motion  pictures  were  made  of  the 
flowfleld,  with  the  f ield-of-view  extending  from  the  nozzle  exit  to  x  -  14.6. 
Spark  photographs  were  taken  with  a  10  ns  spark  source. 
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4.  TEST  PARAMETERS 

In  the  absence  of  excitation,  the  model  displayed  a  one-parameter  family 
of  flows,  obtained  by  adjusting  the  two  exit  doors  to  control  the  exhaust 
area  (A^).  The  parameter  characterizing  members  of  this  family  was  chosen  to 
satisfy  the  unusual  requirements  of  this  experiment. 

The  conventional  parameter  choice  is  the  time-mean  exhaust  area  since  it 
is  directly  related  to  the  inlet  mass  flow  and  varies  monotonically  throughout 
the  available  range  of  flow  conditions.  Unfortunately,  the  present,  unusual 
throttling  arrangement  could  not  be  used  reliably  to  measure  the  mass  flow, 
for  reasons  explained  in  Section  3.2. 

The  ratio  of  the  time-mean  freestream  total  and  exit  static  pressures 
(ve>  is  a  conceptually  sound  parameter  that  could  be  precisely  measured  and 
has  been  used  successfully  in  earlier  work  on  supercritical  diffuser 
flows.  *  However,  in  the  present  inlet  configuration,  vg  becomes  nonunique 
and  therefore  of  limited  value.'  When  the  shock  is  located  at  the  cowl  lip, 
vg  reaches  a  minimum.  Values  greater  than  the  minimum  are  ambiguous;  they 
may  designate  conditions  with  the  mean  shock  position  either  before  or  after 
the  cowl  lip.  Since  this  study  is  focused  on  phenomena  centered  around  just 
this  range  of  ambiguity,  was  not  appropriate  for  the  present  work. 

The  exit  station  Mach  number  (Me)  had  the  fewest  disadvantages  and  was 
adopted  as  the  principal  parameter.  According  to  one-dimensional  isentroplc 
gas  dynamics,  Me  varies  monotonically  with  A  and  is  thus  free  of  the  ambi¬ 


guity  problem  associated  with 
relation  is  double-valued. 


Me  uniquely  defines  ve,  but  the  inverse 


Since  the  flow  at  the  exit  station  was  nonuniform,  an  average  Mach  number 
(Me>  was  defined  from  midspan  profiles  of  time-mean  flow  properties  as 


(a  ,2 . ..!}[ jiv . 

e  yR  /  Pu  Tdy 


The  use  of  this  definition  is  justified  by  its  close  connection  to  the 
Integral  energy  balance  for  the  model.  If  an  integral  energy  equation  is 
written  for  a  control  volume  enclosed  by  the  model  and  bounded  by  a  plane  sur- 


A  •*.  -\ 
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face  at  the  exit  station,  then  the  surface  integral  associated  with  the  exit 
flow  contains  two  terns,  representing  the  kinetic  and  theraal  energy  fluxes. 
M^,  as  defined  by  Equation  (2),  is  proportional  to  the  ratio  of  those  two 
terns.  (The  factor  of  proportionality  is  2/(y  -  1).)  This  definition  was 
preferred  to  a  simple  average  of  Me(y),  which  cannot  be  directly  associated 
with  integral  considerations. 

Local  Mach  numbers  at  the  exit  station  were  determined  from  time-mean 
total  pressure  data  from  a  fixed,  12-tube  rake  located  at  midspan  and  from  the 
time-mean  wall  static  pressure  obtained  as  the  average  of  four  orifices 
located  at  the  same  station. 

was  based  on  midspan  data  only  and  is  greater  than  the  area-average 
Mach  number  by  varying  amounts  because  of  sidewall  boundary-layer  presence. 
This  discrepancy  was  accepted  because  the  measurement  of  Mach  number  distribu¬ 
tions  over  the  entire  exit  area  would  have  been  prohibitively  time-consuming. 
The  discrepancy  is  a  measure  of  deviation  from  two-dimensionality. 

When  excitation  is  used,  the  amplitude  and  frequency  of  the  perturbation 
must  also  be  characterised. 

The  excitation  frequency  (fex)  was  determined  by  the  exit-area  modula¬ 
tion,  occurring  at  twice  the  rotor  shaft  frequency.  The  output  from  the 
shaft-mounted  encoder  provided  a  pulse  train  that  was  counted  over  fixed  time- 
intervals  to  yield  a  precise  measurement  of  f#x. 

The  intensity  of  excitation  was  measured  by  the  exit  amplitude  (Pef)» 
defined  as  the  rms  value  of  the  narrow-band-filtered  wall  static  pressure  sig¬ 
nal  at  the  exit-station  top  wall.  The  filter  center  frequency  was  slaved  to 
the  encoder  output  from  the  exciter;  the  band-width  was  10  Hs. 

The  exit  amplitude  could  be  varied  continuously  through  the  variation  of 
exhaust-area  modulation  amplitude  by  positioning  the  exciter  flaps.  Five  flap 
settings  were  selected  to  cover  the  available  range  of  amplitudes.  In  addi¬ 
tion  to  zero  and  maximum  modulation,  three  intermediate  settings  were  selected 
so  that  at  &c  ■  0.425  and  fex  •  350  Hz,  the  respective  p#^  values  were  0.75, 
0.5,  and  0.25  of  the  maximum  value.  If  ft  and  f__  differed  from  this  arbl- 
trary  choice,  then  the  p#j  values  obtained  at  the  five  flap  settings  generally 
did  not  form  such  a  simple  arithmetic  sequence. 
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The  flaps  could  be  reset  to  within  0.05  mm,  although  such  accuracy  was 
not  required  since  the  parameter  of  interest,  ,  was  always  independently 
measured. 


The  pressure  ratio  (ve)  is  an  easily  measured  mean  flow  parameter. 
Knowledge  of  its  dependence  on  the  parameters  chosen  for  this  work  may  be 
useful  in  relating  the  present  results  to  other  investigations  and  is  given  in 
Figure  5.  (For  later  convenience,  the  reciprocal  of  vg  is  used  in  the  fig¬ 
ure.)  The  principal  connection  is  between  vg  and  but  the  excitation  fre¬ 
quency  and  exit  amplitude  also  have  significant  effects.  Figure  5  also  illu¬ 
strates  that  a  given  vg  may  be  associated  with  two  values,  representing  two 
different  flowfields. 
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5.  TERMINOLOGY 


Ambiguities  and  conflicting  usage  of  terms  concerning  Inlet  flowflelds 
are  common  In  the  literature.  To  avoid  misinterpretation,  the  terminology 
used  in  this  report  will  be  defined  In  detail  at  the  outset,  specifically  con¬ 
sidering  the  experimental  fact  that  the  observed  inlet  flows  were  never 
completely  steady.  Existing  conventions  will  be  followed  whenever  possible. 

5. 1  Instantaneous  Flow  States 

Classification  of  the  instantaneous  flowflelds  will  be  based  on  the 
streamwise  location  of  the  terminal  shock  (Figure  6).  The  terminal  shock  is 
shown  as  a  single,  nominally  normal  shock,  possibly  with  bifurcated  termina¬ 
tions  at  the  walls. 
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Supercritical 


The  flow  Is  said  to  be  In  the  supercritical  state  if  the  terminal  shock 
is  downstream  of  the  cowl  lip,  subcritical  if  it  is  between  the  ramp  and  cowl 
lips,  and  noncritical  if  ahead  of  the  ramp  lip.  Singular  states  occur  if  the 
shock  is  exactly  at  the  cowl  lip  (critical)  or  exactly  at  the  ramp  lip  (in¬ 
cipient  noncritical).  The  conventional  usage  considers  noncritical  states  as 
part  of  the  subcritical  range;  the  distinction  Introduced  here  is  justified  by 
significant  qualitative  differences  between  subcritical  and  noncritical  flows. 

In  certain  circumstances  the  supersonic/subsonic  deceleration  in  inlets 
may  occur  not  through  a  simple  normal  shock,  but  across  a  complex  structure  of 
several  normal  shocks  in  series  (shock  train).  In  such  cases,  the  term  "shock 
position"  becomes  ambiguous,  and  the  above  classification  may  need  to  be 
modified  to  accommodate  the  added  complications.  In  the  present  experiments, 
shock  trains  were  not  observed  near  singular  states. 

5.2  Time-Mean  Flow  Conditions 

Since  natural  oscillations  are  invariably  present,  any  time-mean  flow 
condition  (fixed  Mg)  involves  a  range  of  shock  positions.  Figure  7  illus¬ 
trates  the  peak-to-peak  ranges  found  in  the  present  model.  The  plot  shows 
that  some  Mg  values  may  be  associated  with  more  than  one  type  of  criticality, 
e.g.,  Mg  -  0.265  is  associated  with  a  shock  range  that  includes  both  super- 
and  subcritical  states. 

It  is  customary  and  convenient  to  associate  time-mean  parameters  with  a 
specific  criticality.  This  association  will  be  based  on  the  state  assumed  by 
the  terminal  shock  in  its  upstream-most  position.  In  the  above  cited  example 
the  shock  is  subcritical  in  its  upstream-most  position,  and  therefore  the  flow 
condition  defined  by  M  -  0.265  is  called  subcritical.  Figure  7  illustrates 
the  three  time-mean  criticality  ranges  in  terms  of  M£,  in  accordance  with  the 
present  convention. 

One  might  consider  the  time-mean  shock  position  as  a  basis  for  the  defi¬ 
nition  of  time-mean  criticality;  however,  as  will  be  elaborated  later,  impor¬ 
tant  qualitative  changes  of  behavior  tend  to  be  keyed  to  the  upstream-most 
shock  position.  Also,  experimental  determination  of  the  upstream-most  posi¬ 
tion  is  simpler,  especially  with  solid  sidewall  inlet  models  which  offer  no 
optical  access  to  the  flow  interior. 


Figure  7.  Shock  position  extremes  without  excitation  and  ranges  of 
criticality  for  time-mean  flow  conditions. 


5.3  Oscillation  Modes 

The  flow  exhibited  a  variety  of  complex  oscillations  which  may  be  charac¬ 
terised  in  different  ways.  One  simple,  purely  kinematic  classification  is 
based  on  shock  position  ranges  assumed  during  the  motion. 
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Oscillations  in  which  the  Instantaneous  flow  patterns  remain  in  the  same 
class  of  criticality  will  be  called  pure  super-,  sub-,  or  noncritical  oscilla¬ 
tions.  Oscillations  involving  two  adjacent  criticality  ranges  are  also  pos¬ 
sible  and  have  been  observed  (Figure  7,  0.255  <  M  <  0.276).  These  will  be 
called  dual  modes,  with  the  additional  qualifiers  of  super/sub  (or  sub/non), 
if  the  range  needs  to  be  described  more  precisely. 

Oscillations  may  also  encompass  all  three  ranges  of  criticality 
(Figure  7,  <  0.22)  in  which  case  the  oscillation  is  labeled  a  triple  mode. 

The  above  definitions  of  oscillation  modes  do  not  depend  on  the  cause 
of  the  motion  and  are  therefore  applicable  to  both  natural  and  forced 
oscillations . 


5.4  Buzz 

The  term  buzz  has  been  used  in  the  literature  to  describe  many  types  of 
inlet  oscillations;  therefore,  it  conveys  little  specific  information  other 
than  the  oscillations  are  sustained  spontaneously  and  the  amplitudes  vary  from 
moderate  to  large. 

We  use  the  term  buzz  to  describe  oscillations  in  which  periodic  changes 
of  instantaneous  criticality  occur.  This  Includes  the  dual  and  triple  modes. 
The  triple  mode  will  be  alternately  and  equivalently  referred  to  as  full  buzz, 
since  In  this  mode  the  oscillation  amplitudes  were  much  greater  than  in  any  of 
the  other  modes. 

5.5  Unstart 

Unstart  is  defined  as  a  process  of  a  single  transition  into  a  condition 
that  is  noncritical  in  both  the  time-mean  and  instantaneous  sense.  This  defi¬ 
nition  implies  that  the  oscillations  after  the  transition  (if  any)  are  purely 
noncritical.  The  inlet  remains  in  this  condition  until  some  externally  in¬ 
duced  change  occurs  in  the  boundary  conditions.  The  initial  flow  condition 
and  the  initial  oscillation  may  be  of  any  type  (other  than  noncritical). 

In  contrast  to  buzz,  whose  definition  is  based  on  changes  of  instantane¬ 
ous  flow  states,  unstart  is  defined  primarily  as  a  change  of  time-mean  condi- 


1 


6 


tions.  Full  buzz,  or  dual  (sub/non)  buzz  is  thus  not  to  be  viewed  as  a  peri¬ 
odic  succession  of  unstarts  since  no  sustained  changes  of  time-mean  flow 
conditions  are  involved. 


5.6  Criticality  Boundaries 

The  terminology  established  in  this  section  necessitates  a  reinterpreta¬ 
tion  of  Task  III  outlined  in  Section  2.  The  task  statement  describes  the 
boundary  between  super-  and  subcritical  flow  conditions  and  refers  to  it  as 
unstart  boundary.  This  usage  is  inconsistent  with  Section  5.5  and  will  be 
remedied  as  follows : 

The  set  of  flow  conditions  intermediate  between  two  adjacent  conditions 
of  time-mean  criticality  will  be  referred  to  as  criticality  boundaries.  There 
are  two  such  boundaries:  the  critical  boundary  (CB)  corresponding  to  the 
limiting  case  between  the  super-  and  subcritical  conditions,  and  the  incipient 
noncritical  boundary  (1NB)  separating  the  subcritical  and  noncritical  condi¬ 
tions.  The  boundaries  are  defined  in  terms  of  the  three  time-mean  parameters 

a 

(M  ,  f  ,  and  p  ,)  and  are  themselves  time-mean  entities, 
e’  ex’  *ef' 

Because  of  the  definitions  introduced  in  the  preceding  sections,  Task  III 
is  reinterpreted  to  imply  the  determination  of  the  above  two  criticality 
boundaries . 

The  term  unstart  boundary  (UB)  is  logically  reserved  to  designate  the  set 
of  flow  conditions  causing  unstart  in  the  sense  of  Section  5.5.  The  deter¬ 
mination  of  UB  was  not  possible  in  the  present  experiment  because  noncritical 
shock  patterns  were  dominated  by  the  nozzle  exit  configuration  and  did  not 
resemble  the  detached  shock  pattern  expected  to  exist  in  an  unconfined  free- 
stream  (Figure  6).  The  reliable  determination  of  UB  requires  a  wind  tunnel 
test  or  possibly  a  free  jet  experiment  with  a  large  jet  diameter  compared  to 
the  model  size. 


6.  UNEXCITED  FLOWS 


Reference  9  describes  in  detail  the  behavior  of  the  inlet  model  with 
steady  exit  boundary  conditions.  A  brief  summary  of  this  work  is  included 
here  for  completeness. 

The  procedure  for  operating  the  model  was  to  open  the  doors  fully  and 
pressurize  the  plenum  sufficiently  (at  least  520  kPa)  to  cause  the  initially 
formed  shock  system  to  be  swallowed  by  the  inlet.  The  doors  were  subsequently 
adjusted  to  reach  the  desired  operating  condition.  The  terminal  shock  moved 
upstream  as  the  door  was  closed;  it  could  be  located  anywhere  within  the  in¬ 
let,  and  driven  completely  into  the  nozzle  by  closing  the  doors  sufficiently. 

The  present  contract  work  covered  all  three  ranges  of  criticality,  from 
noncritical  conditions  (Mg  =  0.12)  to  moderate  supercriticality  (Me  *  0.42). 

In  terms  of  a  commonly  used  engineering  definition  of  supercriticality, 

[(pec  -  pg)  /  p  ],  the  maximum  supercriticality  was  18%. 

6.1  Time-Mean  Behavior 

For  the  range  of  flow  conditions  investigated,  the  flow  was  symmetric 
with  respect  to  the  geometric  center  plane  and  possessed  a  reasonably  twa- 
dimensional  central  region.  Figure  8  shows  Mach  number  distributions  at  the 
exit  station  (obtained  in  the  study  of  Reference  9),  deduced  from  data 
obtained  with  a  traversable,  12-pronged  total-pressure  rake.  Two- 
dimensionality  is  expected  to  deteriorate  in  the  streamwise  direction,  so  that 
the  illustrated  distributions  represent  the  worst  case  for  each  flowfield. 
There  was  no  clear  evidence  of  flow  separation  in  the  oil  flow  traces  taken  on 
either  the  top  or  bottom  wall  for  Me  <  0.4. 

The  top-  and  bottom-wall  static-pressure  distributions  are  displayed  in 
Figures  9a  and  9b.  For  Figure  9a,  the  flow  is  fully  supercritical,  and  in 
Figure  9b  the  flow  proceeds  from  critical  to  noncritical  conditions.  The 
supercritical,  bottom-wall  distributions  suggest  the  presence  of  shock-induced 
separation  by  abrupt  changes  of  slope  immediately  behind  the  shock  (x  *  3) 
for  Mfi  >  0.3  (compare  with  Figure  9,  Reference  12).  This  observation  is  in 
apparent  conflict  with  the  absence  of  separation  as  indicated  by  oil  flow 
traces.  Since  oil  patterns  represent  time-mean  information,  the  presence  of  a 
separation  bubble  may  have  been  obscured  by  the  large  amplitude  oscillations. 
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Figure  9.  Wall  pressure  distributions;  (a)  supercritical  and  (b)  critical 
and  subcritical. 

especially  If  the  bubble  was  small.  Separation  Is  not  suggested  by  the 
bottom-wall  pressure  distributions  for  subcritical  conditions,  or  for  the  top 
wall  under  any  conditions.  This  result  Is  unexpected  because  the  terminal 
shock  Mach  number  (Mo)  is  considerably  greater  than  1.3,  the  value  thought 
to  represent  an  approximate  threshold  for  incipient  shock-induced  separa¬ 
tion.^"’^  No  explanation  Is  presently  available. 


6.2  Natural  Oscillations 

As  Illustrated  la  Figure  7,  the  Inlet  exhibited  natural  oscillations 
under  all  conditions.  Pure  supercritical  and  subcritical  oscillations  as  well 
as  a  dual  (super/sub)  mode  were  found  to  exist.  When  the  inlet  was  at  the  INB 
boundary,  a  minute  additional  closing  of  the  doors  caused  the  inlet  to  go  into 
full  buzz  (triple  mode),  l.e.,  noncritical  conditions  were  invariably  associ¬ 
ated  with  full  buzz.  Nc  pure  noncritical  and  no  sub/non  dual  modes  were  ob¬ 
served,  probably  because  of  the  proximity  of  the  nozzle  exit,  which  strongly 
influenced  the  shock  pattern  displayed  by  the  inlet  in  noncritical  states. 

Supercritical  oscillations  appear  to  be  strongly  affected  by  shock/ 
boundary-layer  interactions.  The  oscillations  are  random,  although  the  power 
spectral  density  distributions  of  pressures  may  indicate  one  or  two  broad 
peaks. 

Oscillations  in  which  the  shock  moves  ahead  of  the  cowl  lip  (even  if  tem¬ 
porarily)  tend  to  be  strictly  periodic  and  display  a  spectrum  composed  of  a 
set  of  sharp  peaks  extending  up  to  six  significant  harmonics.  The  fundamental 
frequency  was  near  65  Hz.  In  all  these  cases,  the  vortex  sheet  generated  by 
the  bifurcation  point  (formed  by  the  ramp  shock  and  the  terminal  shock)  is  a 
significant  feature  of  the  flow.  The  vortex  sheet  is  ingested  by  the  inlet 
and  develops  into  a  shear  layer  that  broadens  to  nearly  half  of  the  channel 
height  by  x  2  16. 

Further  details  on  natural  oscillations  can  be  found  in  Reference  9. 


7.  EXCITED  FLOWS 


la  a  real  ramjet,  downstream  perturbations  are  Introduced  by  combustion 
Instabilities  in  the  burner  where  the  burning  depends  on  the  spatial  pressure/ 
velocity  distributions  of  the  flow  entering  the  combustor.  The  inlet  and 
burner  oscillations  are  coupled:  any  perturbation  of  the  entering  flow  changes 
the  combustion  process,  which  in  turn  changes  the  downstream  boundary  condi¬ 
tions  for  the  inlet,  thereby  reacting  back  to  the  inlet  flow.  A  theoretical 
description  of  this  process  requires  the  ability  to  predict  the  response  of 
the  inlet  to  downstream  boundary  conditions  and  the  response  of  the  combustor 
to  the  entering  flow  properties.  The  present  program  Isolates  and  focuses  on 
the  first  of  these  two  problems. 

In  the  present  experiments,  the  combustion  instabilities  are  simulated  by 
modulation  of  the  exit  area  independently  of  the  response  evoked  from  the 
inlet.  This  situation  differs  from  the  Interactive  oscillations  in  ramjets 
with  combustion;  nevertheless,  the  two  systems  can  be  compared  meaningfully  if 
interest  is  confined  to  the  inlet  portions.  The  separation  of  the  Inlet  from 
the  combustor  requires  a  somwewhat  arbitrary  definition  of  a  boundary  surface 
separating  the  two.  In  the  present  experiment,  the  downstream  boundary  of  the 
inlet  was  set  arbitrarily  at  xfi  *  21.91  (exit  station).  If  the  inlet  geometry 
and  the  upstream  boundary  conditions  are  given  (as  they  will  be  throughout 
this  report),  then  spectficaton  of  the  time-dependent  exit  pressure  (or.  If 
nonuniform,  the  distribution  of  exit  pressure)  completely  defines  the  entire 
inlet  flow.  The  manner  in  which  the  exit  pressure  variation  is  created  is 
immaterial,  and  a  mechanically  excited  flow  is  Identical  to  a  combustor- 
excited  flow,  provided  the  histories  of  exit  station  conditions  are  identical. 

Operation  of  the  exciter  created  periodic  perturbations  of  all  flow  pro¬ 
perties.  The  perturbations  propagated  upstream  at  the  speed  of  sound  with 
respect  to  the  fluid.  These  perturbations  triggered  complex  responses  from 
all  points  of  the  inlet  that  either  propagated  or  were  convected  downstream. 
The  perturbation  pressure  at  the  exit  station  is  the  sum  of  the  original  per¬ 
turbation  and  the  net  result  of  all  local  responses;  the  two  contributions 
cannot  be  separated  in  a  steady  oscillation.  The  measured  exit  station  per¬ 
turbations  thus  should  not  be  thought  of  as  inputs;  they  represent  the  input 
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plus  Che  output  returned  by  Che  system.  For  this  reason,  p^  Is  deliberately 
called  the  exit  amplitude  In  preference  to  the  possibly  misleading  term 
excitation  amplitude. 

The  time-dependent  part  of  the  measured  flow  quantities  contained  peri¬ 
odic  contributions  Introduced  by  the  exciter  and  random  contributions  from 
turbulent  fluctuations.  Interest  was  focused  on  the  first  harmonic  of  the 
periodic  contribution,  obtained  by  narrow-band  filtering  of  the  respective 
signals  at  the  excitation  frequency. 

The  parameters  of  the  three-dimensional  test-matrix  were  the  exit  Mach 
number  (M  ),  the  excitation  frequency  (fex)  and  the  exit  (rms)  amplitude 
(pef).  This  set  Is  reasonably  independent  of  facility-related  effects, 
uniquely  characterizes  the  flow,  and  is  acceptable  as  input  to  theoretical 
considerations.  The  experimental  control  of  these  parameters,  however,  was 
neither  direct  nor  independent. 

The  excitation  frequency  was  the  only  variable  that  could  be  set  to  any 
desired  value  independently  of  the  other  two  parameters  by  controlling  the 
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rotor  speed.  In  contrast,  p^  depended  strongly  on  the  setting  of  all  three 
controls:  rotor  speed,  doors,  and  flaps.  depended  primarily  on  the  door 
setting,  but  rotor  speed  and  flap  setting  also  had  nonnegligible  influences 
(Figure  5). 

The  ideal  method  of  executing  a  test  series  having  three  parameters  is  to 
compose  it  of  sub-series  in  each  of  which  two  variables  are  held  constant  and 
only  the  third  is  varied.  In  a  three-dimensional  parameter-space,  the  points 
representing  the  test-conditions  would  then  lie  on  straight  lines  parallel  to 
the  axes.  Because  of  the  indirect  control  of  the  test  parameters  in  the  pre¬ 
sent  experiment,  there  was  concurrent  variation  of  at  least  two  independent 
variables,  so  that  the  representative  points  of  a  sub-series  generally  formed 
a  two-dimensional  (and  sometimes  even  three-dimensional)  curve  in  the  para¬ 
meter  space.  Therefore,  interpolation  is  usually  necessary  to  present  the 
data  in  graphic  form. 

Time  and  cost  constraints  placed  severe  limitations  on  the  achievable 
parametric  resolution.  Ten  values  for  each  of  the  three  parameters  represent 
a  relatively  crude  resolution,  yet  imply  a  massive  test  schedule  of  a  thousand 
runs.  Considering  that  17  time-dependent  flow  properties  could  be  measured. 
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the  results  would  be  recorded  as  17  000  time  traces.  The  complete  exploration 
of  all  properties  as  functions  of  all  parameters  over  their  full  ranges  Is 
thus  not  feasible.  Judgements  and  choices  were  made  In  defining  the  test 
schedule,  attempting  to  focus  on  the  most  critical  phenomena. 

The  mode  of  oscillation  occurring  at  a  specific  flow  condition  could  be 
determined  reliably  from  high  speed  films.  The  potentially  possible  oscilla¬ 
tory  modes  (defined  In  Section  S.3)  are  the  same  with  and  without  excitation, 
but  the  actually  occurring  modes  depended  on  exit  amplitude.  Figure  10  illu¬ 
strates  the  shock  position  ranges  observed  during  forced  oscillations  as  func¬ 
tions  of  Mfi,  superimposed  on  a  plot  of  natural  oscillation  ranges.  All  forced 
oscillation  data  In  the  plot  were  taken  with  the  flaps  at  maximum  setting. 

The  amplitudes  are  greater  by  as  much  as  an  order-of-magnitude  when  excitation 
Is  present.  One  consequence  of  the  large  amplitudes  is  that  purely 
subcritical  oscillations  were  not  found;  only  pure  supercritical,  sub/super 
dual,  and  full  buzz  modes  exist. 

Figure  10  illustrates  that  excitation  enlarges  the  shock  position  ranges 
at  a  given  value  of  mostly  by  extending  them  downstream,  while  the 
upstream-most  extreme  moves  upstream  by  only  a  small  amount.  The  time-mean 
shock  position  moves  downstream  by  up  to  two  throat  heights,  i.e.,  the 
presence  of  excitation  has  the  same  effect  on  the  mean  flow  condition  as  would 
an  increase  of  exhaust  area. 

An  alternative  way  of  expressing  this  observation  is  to  state  that  if  the 
time-mean  shock  position  is  kept  constant,  then  the  imposition  of  forced  os¬ 
cillations  reduces  the  associated  M  .  The  unexcited,  time-mean  inlet  charac- 
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terlstics  are  such  that  lower  Mg  is  associated  with  lower  total-pressure 
losses  in  the  subsonic  region,  suggesting  that  forced  oscillations  alleviate 
subsonic  loss  mechanisms.  No  explanation  is  available  at  present  for  this 
apparent  trend. 
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thus  not  feasible.  Judgeaents  and  choices  were  aade  In  defining  the  test 
schedule,  attempting  to  focus  on  the  most  critical  phenomena. 

The  mode  of  oscillation  occurring  at  a  specific  flow  condition  could  be 
determined  reliably  from  high  speed  films.  The  potentially  possible  oscilla¬ 
tory  modes  (defined  In  Section  5.3)  are  the  same  with  and  without  excitation, 
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The  amplitudes  are  greater  by  as  much  as  an  order-of-magnltude  when  excitation 
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upstream-most  extreme  moves  upstream  by  only  a  small  amount.  The  time-mean 
shock  position  moves  downstream  by  up  to  two  throat  heights,  i.e.,  the 
presence  of  excitation  has  the  same  effect  on  the  mean  flow  condition  as  would 
an  increase  of  exhaust  area. 

An  alternative  way  of  expressing  this  observation  is  to  state  that  if  the 
time-mean  shock  position  is  kept  constant,  then  the  imposition  of  forced  os¬ 
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terlstics  are  such  that  lower  is  associated  with  lower  total-pressure 
losses  in  the  subsonic  region,  suggesting  that  forced  oscillations  alleviate 
subsonic  loss  mechanisms.  No  explanation  is  available  at  present  for  this 
apparent  trend. 
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Figure  10.  Shock  position  extremes  for  forced  oscillntions  with 
maximum  flap  setting.  Extremes  for  natural  oscillations  are  also 
shown  for  comparison. 


7.1  Supercritical  Oscillations 

In  this  class  of  motions  the  shock  always  stayed  downstream  of  the  cowl, 
the  flow  at  the  geometric  throat  was  always  steady,  and  the  mass  flow  into  the 
inlet  was  constant.  These  conditions  combined  to  create  a  relatively  simple 
type  of  motion. 

Figure  11  illustrates  how  the  exit  amplitude  varies  with  excitation  fre¬ 
quency  if  the  flap  setting  is  held  constant.  The  concurrent  variation  of  the 
two  principal  parameters  is  evident.  Results  from  (filtered)  wall  pressure 
transducers  located  at  interior  points  yield  similar  plots,  but  with  major 
differences  in  the  locations  and  magnitudes  of  the  minimum  or  maximum  values. 
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In  light  of  the  concurrent  variation  of  fgx  and  pgf,  it  is  difficult  to 
separate  the  frequency  and  amplitude  effects. 
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Figure  11.  Dependence  of  exit  amplitude  on  excitation  frequency 
for  four  different  flap  settings.  Me= 0.420. 


The  difficulties  were  greatly  alleviated  by  the  fact  that  the  internal 
(wall)  pressure  perturbation  amplitudes  varied  linearly  with  the  exit  ampli¬ 
tude.  This  linearity  is  illustrated  in  Figure  12.  Data  taken  at  the  same 
excitation  frequency  fall  closely  on  a  straight  line:  the  correlation  coeffi¬ 
cients  for  least-squares  fit  straight  lines  were  typically  over  0.99.  Linear 
behavior  was  found  for  all  supercritical  conditions,  for  all  frequencies/ 
amplitudes,  and  at  all  locations  where  the  flow  was  subsonic  at  all  times. 
Linearity  failed  only  at  those  locations  that  were  Intermittently  passed  over 
by  the  shock  during  its  downstream  excursions. 

The  slope  of  the  linear  relations  (P)  is  actually  the  local  wall  pressure 
amplitude  normalized  by  the  exit  pressure  amplitude.  P  depends  on  Mg  and  the 
excitation  frequency  but  not  on  the  exit  amplitude,  which  was  effectively  eli¬ 
minated  as  a  variable  by  normalization.  The  number  of  independent  variables 
for  this  class  of  oscillations  is  thus  reduced  from  three  to  two. 

Figure  13  contains  curves  of  normalized  wall  pressure  as  a  function  of 
frequency  for  transducer  T3  (x  -  5.76),  which  is  downstream  of  the  shock. 

Data  for  four  Me  values  are  given.  The  symbols  are  test  results,  while  the 


Figure  12.  Linearity  of  induced  wall  pressure 
fluctuation  amplitudes  for  Me»  0.410  and  *  =  5.76. 
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Figure  13.  Normalized  wail  pressure  fluctuation  amplitudes 
at  x  ■  5.76  for  four  flow  conditions. 


lines  are  somewhat  conjectural  indications  of  what  the  behavior  is  thought  to 
be  between  the  test  points.  The  curves  indicate  a  strong  frequency  depen¬ 
dence,  even  if  conservative  error  bands  are  superimposed. 

The  curves  of  Figure  13  were  used  to  construct  the  contour  plot  of  Figure 
14,  which  includes  additional  interpolations  with  respect  to  Me»  Similar 
contour  plots  were  prepared  for  three  other  transducer  locations  (x  ■  7.37, 
8.99,  and  13.83),  all  of  which  were  found  to  have  nearly  identical  structures. 
It  is  therefore  probable  that  the  contour  plot  of  Figure  14  reflects  actual 
trends,  despite  the  smoothing  and  interpolation  employed.  The  following 
trends  are  recognized: 

•  At  a  near-critical  condition  (Mg  *  0.27),  the  normalized  wall  pressure 
amplitude  peaks  at  approximately  280  Hz.  This  peak  frequency 
gradually  decreases  to  150  Hz  as  is  increased  to  0.33. 

•  As  is  increased  above  0.33,  two  well-defined  peaks  emerge,  one  at 
90  Hz  and  another  at  190  Hz.  These  peak  frequencies  stay 
approximately  constant  up  to  the  experimental  limit  of  Mg  =  0.41. 

To  the  best  of  the  authors'  knowledge,  the  present  results  constitute  the 
first  experimental  demonstration  of  linearity  of  the  response  of  such  a  com¬ 
plex  flowfield  to  relatively  large  perturbations.  Acoustic  theories  assume 
such  linearity,  but  the  amplitude  limitation  on  the  assumption  is  always  a 
question.  In  the  present  model,  linear  response  was  found  at  exit  amplitudes 
up  to  7.5%  of  the  mean  exit  static  pressure.  The  largest  observed  internal 
pressure  amplitude  was  approximately  16%  of  the  exit  static  pressure  (occurred 
at  conditions  other  than  those  causing  the  previously  mentioned  7.5%  exit 
amplitude). 

Demonstration  of  the  linear  system  response  and  determination  of  flow 
conditions  associated  with  peak  sensitivity  were  made  possible  by  the 
variable-amplitude  feature  of  the  exciter  mechanism,  fully  justifying  the 
considerable  effort  expended  in  its  development. 

7.2  Buzz  with  Excitation 

Forced  oscillations  at  sub-  or  noncritlcal  mean  flow  conditions  invari¬ 
ably  Involved  periodic  changes  of  criticality  states  and  are  therefore  classi¬ 
fied  as  buzz.  Dual  and  triple  buzz  both  occurred  (Figure  10). 
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Figure  14.  Contour  plot  of  normalized  pressure  fluctuation  amplitudes  (P)  at  x  =  5.85,  contours  at 
increments  of  0.5.  Dashed  lines  indicate  frequencies  for  maximum  response. 


Triple  buzz  displayed  greater  shock  motion  amplitudes  than  dual  buzz,  but 
the  differences  between  the  two  types  were  not  as  marked  as  they  were  without 
excitation.  Excitation  moved  the  mean  shock  position  downstream,  on  occasion 
changing  a  triple  mode  into  a  dual  one  by  eliminating  the  noncritical  portion 
of  the  oscillatory  range.  The  similarities  allow  joint  discussion  of  the  two 
modes  of  forced  buzz. 


The  case  of  M  =  0.180  with  fixed  maximum  flap  setting  was  selected  for 
detailed  study.  At  this  condition  the  unexcited  flow  displayed  full  buzz  at 
64  Hz.  High-speed  movies  and  top-wall  static-pressure  fluctuation  power  spec¬ 
tral  density  distributions  (PSD's)  were  obtained  without  and  with  excitation 
at  several  frequencies.  Attention  was  focused  on  transducer  T2,  located  at 
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x  =  4.15,  a  location  sufficiently  downstream  of  the  shock  to  remain  subsonic 
at  all  times,  yet  sufficiently  close  to  the  shock  to  assure  similarity  between 
pressure  and  shock  displacement  spectra  at  low  frequencies. 

The  spectrum  and  waveform  of  the  exit  pressure  in  natural  buzz,  pe(t),  is 
shown  in  Figure  15a,  b.  The  fluctuation  is  clearly  dominated  by  periodic  con¬ 
tributions,  consisting  of  several  significant  harmonics  which  are  also  evident 
from  the  distorted  waveform.  Random  contributions  are  negligible. 
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Figure  15.  Pressure  signals  from  top  wall  pressure  transducer  at  x  =  4.15.(a),  (b) 
Me  =  0.180,  noncritical,  unexcited,  full  buzz;  (c)  Me  =  0.430,  supercritical,  excited 
at  148  Hz. 
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The  perturbation  introduced  by  the  exciter  cannot  be  isolated  either  from 
natural  oscillations  or  from  the  response  generated  by  the  inlet.  As  an  ap¬ 
proximation,  the  spectrum  of  pg(t)  for  a  supercritical  flow  (Mg  -  0.430) 
is  shown  in  Figure  15c.  Under  this  condition,  reflections  from  the  terminal 
shock  are  presumably  weak  (Reference  13),  and  the  upstream-moving  wave  intro¬ 
duced  by  the  exciter  is  probably  the  principal  contributor  to  the  measured 
Pe(t).  The  spectrum  is  dominated  by  a  sharp  peak  with  small  higher  harmonics, 
i.e.,  the  excitation  is  essentially  sinusoidal. 

In  this  class  of  motions,  the  natural  and  forced  contributions  do  not 
merely  coexist,  but  interact  with  each  other  in  various  complex  ways.  Figure 
16  shows  a  sequence  of  spectra  obtained  at  4  excitation  frequencies  from  28  to 
86  Hz,  a  range  containing  the  natural  buzz  frequency  (f.  j  =  64  Hz).  The  spec- 
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Figure  16.  Pressure  fluctuation  PSD’s  at  Mt  =  0.180  from  top-wail  transducer  at 
x-4.15,  for  frequencies  near  the  62  Hz  natural  buzz  frequency  (fb,)-  Note  the  presence 
of  numerous  peaks  if  excitation  frequency  (fn)  is  far  from  fb,.  If  f(X  is  close  to  fb,, 
a  simple  set  of  harmonics  is  observed. 


tra  for  f  =  38  and  76  Hz  are  similar  to  the  natural  buzz  spectra  (Figure 
15a);  the  only  difference  is  that  the  fundamental  peak  occurs  at  the  excita¬ 
tion  frequency.  High-speed  films  show  little  difference  between  natural  and 
excited  flows  in  these  two  cases.  The  spectra  for  28  and  86  Hz,  however, 
display  a  significantly  different  structure,  consisting  of  many  more  peaks 
than  either  of  the  two  basic  spectra  (Figure  15a  and  15c). 

High-speed  films  of  forced  buzz  in  these  cases  clearly  show  the  presence 
of  many  contributory  frequencies  and  the  resulting  complex  wave  shapes.  Since 
some  of  the  frequencies  may  be  low,  the  periods  are  long,  not  clearly  recog¬ 
nizable,  and  the  motion  has  the  appearance  of  being  random. 

Further  investigation  of  these  differences  revealed  that  two  distinct 
types  of  interactions  may  occur  between  natural  and  forced  motions;  these  are 
discussed  in  the  next  two  sections. 

7.2.1  Modification  of  Buzz 

If  the  excitation  frequency  was  close  to  the  natural  buzz  frequency 
(fbl),  then  the  resultant  motion  was  similar  to  unexcited  buzz;  it  merely 
occurred  at  fex,  at  larger  amplitudes.  The  shock  position  range  was  displaced 
downstream,  which  in  some  cases  changed  the  mode  of  oscillation  from  triple 
buzz  to  dual  buzz.  This  modification  of  the  natural  mode  occurred  over  the 
frequency  range 

f 

0.5  <  <  1.3  (3) 

rs.  T  />. 

bl 

and  is  illustrated  by  the  38  and  76  Hz  cases  in  Figure  16. 

Modification  also  occurred  when  fgx  was  close  to  a  higher  harmonic  of 
f^j,  although  the  frequency  range  for  such  coupling  was  increasingly  narrower 
for  higher  harmonics.  In  such  cases  the  fundamental  frequency  was  a  submul¬ 
tiple  of  the  excitation  frequency 

f,  ■  f  /  n,  n  *  2,  3...  . 

1  ex 
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fj  was  always  close  to,  but  not  necessarily  equal  to,  fbj.  Figure  17 
Illustrates  such  forced  buzz  spectra  for  the  n  *  2,  4,  5,  and  6  cases. 

If  modification  occurs  at  the  nfc**  submultiple  of  fex,  then  the  amplitudes 
of  the  n^  peak  and  of  the  fundamental  both  increase,  while  the  other  peaks 
remain  largely  unaffected.  Furthermore,  the  greatest  amplitude  always  occurs 
at  the  fundamental,  regardless  of  which  harmonic  couples  to  the  excitation 
frequency.  This  observation  indicates  a  transfer  of  energy  among  harmonics, 
which  is  a  typically  nonlinear  phenomenon. 

7.2.2  Modulation 

If  the  excitation  frequency  is  not  close  to  the  natural  buzz  frequency  or 
one  of  its  harmonics,  then  a  different  type  of  interaction  takes  place:  the 
excitation  modulates  the  buzz  (for  fex  <  fbj)  or  the  buzz  modulates  the 
excitation  (for  fbl  <  fex)* 


Figure  18  shows  the  PSD  for  the  case  where  fex  is  between  and  f ^ • 

The  natural  buzz  pattern  (fbn  -  nfbl,  n  -  1,  2,  3..*)  appears  with  the  peak 
representing  the  excitation  (f  »  f  ) .  In  addition,  a  number  of  new  peaks 
appear,  whose  frequencies  are  found  to  be  given  precisely  by 

fj  “  |fex  *  jfbl‘»  j  "  lt  2*  3***  »  (4) 

i.e.,  the  frequencies  are  given  as  the  absolute  values  of  the  sums  and  differ¬ 
ences  of  the  excitation  frequency  and  the  buzz  harmonics. 

The  origin  of  these  frequencies  can  be  found  by  considering  the  trigono¬ 
metric  identity, 

2  cos  (2wfext)  cos(2*jfblt)  -  cos[2x(fex  - 

+  cos[2it(fex  +  jfbl)t],  (5) 

from  which  it  is  evident  that  the  measured  waveform  can  be  described  as  a 
product  of  two  basic  waves.  One  basic  wave  is  related  to  the  excitation, 

Pex(t)  -  1  +  Ex  cos(2nfext)  ,  (6) 

and  the  other  describes  the  buzz  and  its  harmonics. 
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fjgare  It.  Modulation  of  buzz  by  excitation. 

M#«  0.180,  x*4.1S. 


(7) 


6 

Pb(t)  -  Bj  cos(2njfblt), 

j-1 

where  Ej  and  Bj  are  constant  coefficients. 

The  process  of  generating  a  new  wave  by  multiplying  two  basic  waves  is 
called  modulation  of  the  higher-frequency  wave  by  the  lower-frequency  one.  In 
some  cases,  sidebands  occur  on  both  sides  of  the  excitation  peaks  (Figure  19) 
that  are  reminiscent  of  spectra  discussed  in  amplitude-modulated  radio  trans¬ 
mission  technology. 

Examples  of  modified  and  modulated  pressure  signals  are  compared  in  Fig¬ 
ure  20.  The  modified  wave  shows  the  evidence  of  higher  harmonics,  but  it  is 
dominated  by  the  fundamental.  Successive  waveforms  are  repeated  for  each 
fundamental  period.  The  modulated  wave  shape  is  much  more  complex,  in  part 
due  to  the  fact  that  there  are  buzz  harmonics  both  below  and  above  the 
excitation  frequency,  and  there  is  no  clearly  defined  carrier  wave.  The 
essential  feature  (that  the  measured  wave  is  a  product  of  two  other  waves)  can 
only  be  verified  from  the  PSD's. 
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Figure  19.  Illustiition  of  sidebands  introduced 
by  modulation,  Me- 0.180,  x>4.15. 
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Cjgnre  20.  Modified  ud  modulated  waveforms, 

Mc =0.180,  if ■  4.15.  (a)  f„  =  346  and  (b)  f„  =  150. 


7.2.3  Linearity 

The  discussion  of  the  previous  sections  clearly  Indicates  that  both  types 
of  forced  buzz  are  highly  nonlinear  phenomena.  Linear  amplitude  response, 
which  so  conveniently  simplified  the  description  of  purely  supercritical  os¬ 
cillations,  is  not  expected.  Even  if  it  were  found,  it  would  be  accidental 
since  no  linear  theory  can  describe  the  observed  spectral  transfer  of  energy 
and/or  the  amplitude  modulation  observed. 

Correspondingly,  no  attempt  was  made  to  establish  the  response  vs.  exci¬ 
tation  amplitude  characteristics  in  this  mode. 
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8.  CRITICALITY  BOUNDARIES 

8.1  Procedures 

Passage  of  the  terminal  shock  over  the  cowl  or  the  ramp  lips  was  detected 
by  fast-response  pressure  transducers  located  on  the  outside  of  the  cowl  and 
ramp  surfaces,  i.e.,  on  the  downstream  walls  of  the  large  suction  slots  separ¬ 
ating  the  nozzle  from  the  inlet  (Figure  3).  For  supercritical  conditions,  the 
transducer  signals  contained  low-level  noise  only  (Figure  21a),  while  entry 
into  the  subcrltlcal  state  produced  disturbances  in  the  cowl  slot  signal, 
illustrated  in  Figure  21b.  Entry  into  the  noncrltical  state  was  invariably 
associated  with  full  buzz  and  the  abrupt  appearance  of  square-wave  like  sig¬ 
nals  from  the  ramp  slot  transducer  (Figure  21d). 

The  boundaries  were  approached  by  gradually  closing  the  doors  while 
keeping  excitation  frequency  and  flap  position  constant.  The  slot  transducer 


t  (ms)  t  (ms) 

GPJI4II4-M 

Figure  21.  Pressure  traces  from  transducers  located  in  ramp  slot  (RS)  and  cowl  slot  (CS). 
(a)  Supercritical,  (b)  critical,  (c)  subcrltlcal,  and  (d)  noncrltical  (full  buzz). 
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signals  were  monitored  visually  during  the  process;  the  first  appearance  of 
perturbations  indicated  that  the  boundary  condition  had  been  reached.  The 
door-closing  process  was  associated  with  variations  of  the  excitation  ampli- 

A  A 

tude,  Peft  so  that  the  value  of  pe^  associated  with  the  boundary  was  a  test 
result  rather  than  a  controlled  variable. 

Figures  22,  23  and  24  illustrate  the  boundaries.  Figure  22  shows  the 
limiting  values  as  a  function  of  frequency  for  four  different  flap  posi¬ 
tions,  and  Figures  23  and  24  show  the  pe^  values  measured  at  the  boundaries  in 
the  same  test  series.  All  symbols  are  experimental  data.  Data  taken  with  the 
same  flap  settings  are  connected  by  straight  lines  for  easier  visualization. 


Figare  22.  M(ui  reaction  of  f„  for  various  flap  lettings  at  the  critical  and  the  Incipient 
aoacrHkal  boundaries. 


Figures  22,  23  and  24  represent  the  experimental  Information  available  in 
a  fora  closely  related  to  the  experimental  procedures.  For  interpreting  the 
results,  the  flap  setting  is  of  no  importance  and  ms  eliminated  by  combining 
the  two  figures  to  produce  plots  of  Mg  values  as  a  function  of  excitation  in¬ 
tensity  and  frequency  (Figures  25  and  26).  These  figures  contain  only  opera¬ 
tional  parameters  of  direct  concern,  at  the  cost  of  somewhat  reuuced  precision 
introduced  by  the  interpolations  needed  in  the  cross-plotting  procedure. 

8.2  Discussion  of  Criticality  Boundaries 

Figures  25  and  26  represent  the  boundaries  referred  to  in  Section  2 
(Task.  Ill)  and  more  precisely  defined  in  Section  5.6.  They  represent  perhaps 
the  most  useful  engineering  information  in  this  report. 

A 

The  shape  of  the  explored  (p  ,,  f  )  domain  is  irregular  in  each  case, 

A  W  W 

because  the  maximum  pg^  value  obtainable  by  keeping  the  flaps  at  maximum 
setting  was  a  strong  function  of  ffiX.  The  domain  shape  thus  reflects  experi- 


Figure  23.  p^  as  a  function  of  fH  for  various  flap  settings  at  the  critical  boundary. 
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mental  limitations  only  and  is  therefore  largely  irrelevant.  Given  time  and 
resources,  the  boundaries  could  be  determined,  in  principle,  for  additional 

A 

regions  of  the  Pef(fex)  plane. 

The  important  information  is  contained  in  the  values  of  which  is  pro¬ 
portional  to  the  degree  of  supercriticality.  Higher  Mg  Indicates  that  a 
higher  margin  is  required  to  prevent  reaching  the  critical  boundary  (CB)  or 
the  incipient  noncritical  boundary  (INB).  Since  the  1NB  is  also  the  boundary 
for  the  onset  of  full  (triple)  buzz,  INB  is  probably  the  more  important  of  the 
two. 


A  dominant  feature  shared  by  the  two  boundaries  is  the  rapid  increase  of 
Me  as  frequencies  decrease  below  approximately  75  Hz.  This  frequency  cor¬ 
responds  to  the  upper  limit  at  which  the  excitation  can  couple  to  the  buzz 


fundamental  (Section  7.2.1).  Frequency  dependence  is  weak  above  75  Hz. 


Dependence  on  exit  amplitude  is  different  for  the  two  boundaries.  For 

the  CB,  the  M  value  at  the  boundary  peaks  at  p  ,  -  2  5  kPa  regardless  of 
e  er 

frequency.  This  feature  is  puzzling  since  the  expectation  is  that  M  should 

A  C 

increase  monotonically  with  pg^:  increasing  excitation  amplitudes  will  be 
associated  with  increasingly  greater  shock  displacements,  and,  since  the 
upstream-most  position  (minimum  x^)  is  held  fixed  at  the  cowl  or  the  ramp,  the 
mean  shock  position  will  move  downstream.  In  the  absence  of  excitation,  such 
a  position  is  associated  with  stronger  shocks,  higher  losses,  and  a  higher 

A 

exit  Mach  number.  The  experimental  trend  for  pgf  >  5  kPa  is  the  opposite  of 
this  expectation.  One  possible  cause  of  this  anomaly  might  be  the  reduction 
of  subsonic  viscous  losses  with  increasing  oscillation  amplitude,  as  already 
mentioned  in  Section  7. 

The  amplitude  dependence  of  the  INB  (Figure  25)  differs  considerably  from 
that  found  in  the  CB.  The  INB  displays  little  amplitude  dependence  above  the 
buzz  frequency,  but  below  it,  3  increases  sharply  and  linearly  with  p  . 

In  combination  with  the  already  described  frequency  dependence,  the  low- 
frequency/  high-amplitude  combination  is  clearly  the  most  likely  to  cause  the 
inlet  to  reach  the  INB  and  thereby  the  onset  of  full  buzz. 

The  subcritical  conditions  immediately  adjacent  to  the  INB  correspond 
mostly  to  dual  buzz  mode  oscillations,  with  the  exception  of  low  amplitudes 
and  certain  narrow  ranges  in  which  pure  subcritical  oscillations  may  occur. 
As  explained  in  Section  7.2,  the  excitation  may  either  modulate  or  modify  the 
natural  buzz.  The  ranges  for  each  type  of  coupling  have  not  been  explored  in 
sufficient  detail  to  allow  their  graphic  delineation  to  be  superimposed  on 
Figure  26.  However,  the  absence  of  any  kind  of  anomaly  from  the  boundaries  at 
any  of  the  higher  buzz  harmonics  shows  that  excitation  at  some  multiple  of  the 
buzz  frequency  is  ineffective  and  does  not  demand  an  increase  of  margin  to 
maintain  stable  operation. 
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9.  SUMMARY 


The  response  of  an  external-compression,  M^  =  1.84  inlet  to  periodic 
downstream  perturbations  was  explored  experimentally  over  all  ranges  of  criti¬ 
cality.  The  peak-to-peak  amplitude  of  the  exit-station  pressure  fluctuations 
ranged  up  to  8%  of  the  mean  exit  pressure  at  frequencies  from  20  to  360  Hz. 

The  findings  are  described  using  a  set  of  terms  defined  to  be  valid  in  the 
presence  of  natural  or  externally  induced  unsteadiness. 

In  pure  supercritical  oscillations,  the  natural  fluctuations  are  broad¬ 
band  and  do  not  appear  to  couple  strongly  to  the  enforced  periodic  motion. 

The  amplitudes  of  internal  pressure  fluctuations  in  the  subsonic  region  (fil¬ 
tered  at  the  excitation  frequency)  vary  linearly  with  the  exct  amplitude  under 
all  supercritical  conditions.  The  response  depends  strongly  on  frequency  and 
on  mean  flow  conditions,  and  distinct  frequencies  of  maximum  sensitivity  were 
found. 

Subcritical  and  noncritical  conditions  are  associated  with  natural  oscil¬ 
lations  at  well-defined  discrete  frequencies,  with  minor  contributions  from 
random  (turbulent)  fluctuations.  Downstream  perturbation  at  frequencies  close 
to  the  natural  mode  will  cause  the  natural  mode  to  occur  at  the  excitation 
frequency.  If  the  natural  and  forcing  frequencies  are  too  far  apart,  then  the 
observed  pressure  vs.  time  waveforms  are  precisely  described  as  the  product  of 
the  excitation  waveform  and  the  natural  oscillation  waveform,  including  higher 
harmonics  of  the  latter,  up  to  the  sixth.  The  two  oscillations  modulate  each 
other. 

Criticality  boundaries,  defined  as  sets  of  flow  conditions  separating  the 
supercritical/subcritical  and  the  subcritical/noncritical  flow  condition 
ranges,  were  determined.  The  time-mean  exit  Mach  number  describing  each 
boundary  was  independent  of  frequency  over  75  Hz  and  displayed  sharp  depend¬ 
ences  on  both  amplitude  and  frequency  below  75  Hz. 

The  investigated  flowfields  form  a  heterogeneous  group,  with  large  varia¬ 
tions  of  principal  features  and  dominant  mechanisms.  It  is  clear  that  a  more 
thorough  investigation  is  needed  focusing  on  a  more  narrowly  defined  class  of 
flows,  permitting  a  better  characterization  of  the  physics.  The  principal  ac¬ 
complishment  of  the  present  study  is  to  provide  a  basis  for  the  judicious 
selection  of  meaningful  and  manageable  problems  for  future  inlet  flow 
investigations. 
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